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Empirical tight binding(ETB) methods are widely used in atomistic device simulations. Tradi¬ 
tional ways of generating the ETB parameters rely on direct fitting to bulk experiments or theoretical 
electronic bands. However, ETB calculations based on existing parameters lead to unphysical re¬ 
sults in ultra small structures like the As terminated GaAs ultra thin bodies(UTBs). In this work, 
it is shown that more reliable parameterizations can be obtained by a process of mapping ab-initio 
bands and wave functions to tight binding models. This process enables the calibration of not only 
the ETB energy bands but also the ETB wave functions with corresponding ab-initio calculations. 

Based on the mapping process, ETB model of Si and GaAs are parameterized with respect to hy¬ 
brid functional calculations. Highly localized ETB basis functions are obtained. Both the ETB 
energy bands and wave functions with subatomic resolution of UTBs show good agreement with 
the corresponding hybrid functional calculations. The ETB methods can then be used to explain 
realistically extended devices in non-equilibrium that can not be tackled with ab-initio methods. 

PACS numbers: 


I. INTRODUCTION 

Modern semiconductor nanodevices have reached criti¬ 
cal device dimensions in the sub-10 nanometer range. 
These devices consist of complicated two or three di¬ 
mensional geometries and are composed of multiple ma¬ 
terials. Confined geometries such as ultra thin body 
(UTB)PF inFETs 2 and nanowire^ structures are usually 
adopted in nanometer scale device designs to obtain de¬ 
sired performance characteristics.Most of the electrically 
conducting devices are not arranged in infinite periodic 
arrays, but are of finite extent with contacts controlling 
the current injections and potential modulation. Typi¬ 
cally, there are about 10000 to 10 million atoms in the 
active device region with contacts controlling the current 
injection. These finite sized structures suggest an atom¬ 
istic, local and orbital-based electronic structure repre¬ 
sentation for device level simulation. Quantitative de¬ 
vice design requires the reliable prediction of the mate¬ 
rials’ band gaps and band offsets within a few meV and 
important effective masses within a few percent in the 
geometrically confined active device regions. ETB model 
is usually fitted to bulk dispersions without any defini¬ 
tion of the spatial wave function details. However, recent 
ab-initio study of UTBi^ showed that the surface carrier 
distribution in confined systems is strongly geometry and 
material dependent. This suggests that the charge dis¬ 
tribution for realistic predictions of nanodevice perfor¬ 
mances should be resolved with subatomic resolution. 

Ab-initio methods offer atomistic representations with 
subatomic resolution for a variety of materials. However, 
accurate ab-initio methods, such as Hybrid functionals^ 
G’W® and BSE approximation^ are in general computa¬ 
tionally too expensive to be applied to systems containing 
millions of atoms. Furthermore, those methods assume 
equilibrium and cannot truly model out-of-equilibrium 


device conditions where e.g. a large voltage might have 
been applied to drive carriers. The ETB methods are 
numerically much more efficient than Ab-initio methods. 
ETB has established itself as the standard state-of-the- 
art basis for realistic device simulations. It has been suc¬ 
cessfully applied to electronic structures of millions of 
atoms 8 ^ as well as on non-equilibrium transport problems 
that even involve inelastic scattering!® The accuracy of 
the ETB methods depend critically on the careful cali¬ 
bration of the empirical parameters. The traditional way 
to determine the ETB parameters is to fit ETB band 
structures to experimental data of bulk materials:- 1 - 0 11 

The ETB basis functions remain implicitly defined dur¬ 
ing traditional fitting processes. The lack of explicit basis 
functions makes it difficult to predict wave function de¬ 
pendent quantities like optical matrix elements with high 
precision. More importantly, ETB models parameterized 
by traditional fitting processes suffer from potential am¬ 
biguity when applied to ultra small structures such as 
UTBs, nanowires and more complicated geometries. For 
instance, the existing ETB parameters of GaA^H: applied 
to a As terminated GaAs UTB with an implicit Hydro¬ 
gen passivation modeP 2 results in unphysical top valence 
band states as shown in Fig. [I] The real space probabil¬ 
ity amplitudes of ab-initio topmost valence bands corre¬ 
spond to confined states with the probability amplitude 
peaking in the center of the UTB rather than the surface 
of the UTB as in ETB.In Fig. [I] the hybrid functional cal¬ 
culations include Hydrogen atoms explicitly whereas the 
ETB calculations include only their impact implicitly. 12 
The mismatch between the envelopes of ETB and ab- 
initio wavefunctions suggests a calibration of wave func¬ 
tions in the ETB parameterization process is necessary. 
It is also found that the method of passivation (i.e. im¬ 
plicit or explicit inclusion of Hydrogen atoms) has an 
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effect on the nature of the valence band states. 
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FIG. 1: In As terminated GaAs UTBs, hybrid functional 
probability amplitudes (a) of the top valence bands are con¬ 
fined states with probability amplitude peaking in the center 
of the UTB, while the ETB valence states (b) are surface 
states. 

Therefore, a more fundamental fitting process that 
relates both the band structure and the wave func¬ 
tions of ETB models with ab-initio calculations is 
desirable. Existing approaches to construct local¬ 
ized basis functions and tightbinding-like Hamiltonians 
from ab-initio results inclu de maximally localized Wan- 
nier functions(MLWF) J121M quasi-atomic orbitalsp^ 6 or 
DFT-TB analysis P^The MLWFs are constructed using 
Bloch states of either isolated bands 13 or entangled 
bandsf^These methods typically include interatomic in¬ 
teractions beyond first nearest neighbors. However, these 
methods do not eliminate the above discussed ambiguity 
of the commonly used orthogonal sp3d5s* ETB models. 
Furthermore, these approaches usually disregard excited 
orbitals which are often needed to correctly parameterize 
semiconductors conduction bands. In previous work, it 
was already suggested how to generate ETB parameters 
that are compatible with typical ETB models and still re¬ 
produce ab-initio results.^ 8 This previous method was al¬ 
ready applied to several materials such as GaAs, MgdMl 
and SmSe® and yielded a good agreement between bulk 
ETB and ab-initio band structures. However, the result¬ 
ing wave functions did not satisfactorily agree with the 
ab-initio wave functions. 

In this paper, an improved algorithm of Ref ll8l is pre¬ 
sented that ’’maps” ab-initio results (i.e. eigenenergies 
and eigenfunctions) to tight binding models. Compared 
with the previous workpS the presented method allows 
much better agreement of the ETB and ab-initio wave 
functions. In this present mapping algorithm, rigorous, 
wavefunction-derived ETB parameters for the Hamilto¬ 
nian, for highly localized basis functions, and for ex¬ 
plicit surface passivation are obtained. It is important 
to mention that the ETB Hamiltonian of this method 
can be limited to first nearest neighbor interaction. The 
mapping process is applied to both bulk Si and GaAs 
to generate ETB parameters and explicit basis functions 
from corresponding hybrid functional calculations. It is 
demonstrated in this work, that the wave-function de¬ 
rived ETB Hamiltonian does not yield the ambiguity dis¬ 


cussed with Fig. [T| In the same way, the transferability 
of the ETB model to nanostructures is improved. This 
is demonstrated by a comparison of ETB and Hybrid 
functional results in GaAs and Si UTBs. 

This paper is organized as follows. In section [Hj the 
algorithm of parameter mapping from ab-initio calcula¬ 
tions to tight binding models is described. Section m 
shows the application of the mapping algorithm to bulk 
and UTB systems. Subsection |III A presents the applica¬ 
tion of the present algorithm to bulk Si and GaAs. Bulk 
band structures and realspace basis functions are shown 
and discussed there as well. Subsection IIII Bl shows the 
application of the algorithm to UTB systems and com¬ 
pares ETB band structures and wave functions with cor¬ 
responding ab-initio results. The algorithm and its re¬ 
sults are summarized in Section IIVI 


II. METHOD 

A. Parameter Mapping Algorithm 

The algorithm of the parameter mapping from ab-initio 
results to ETB models is shown in Fig. [2] As will be shown 
in the following, the ETB parameters and basis functions 
are obtained in an iterative fitting procedure that spans 
over 5 steps (with steps 3 through 4 being iterated). The 
resulting 1st nearest neighbor Hamiltonian H TB (k) is of 
Slater Roster table typePSUB The resulting basis ©final 
is composed of orthonormal real space functions ©final = 
{’I'nTm (r) j which have the shape (vectors are given in 
bold type) 

(**) — ($, 0) ^a,n,l (^) “b 

Y. Yi',m' ( 0 , 4>) Ra,n,l,V,m' (r) .( 1 ) 
l ,m 
(l' 

Here, a labels the atom type, whereas the n, l and m 
are principle, angular and magnetic quantum numbers, 
respectively. All materials considered in this work con¬ 
tain no magnetic polarization. Therefore, the basis func¬ 
tions are spin independent. The tesseral spherical har¬ 
monics Yi, m {0, </>) describe the dependence of the basis 
functions on the angular coordinates 6 and </>. The func¬ 
tions R a ,n,i {r) and R a ,n,i,l',m' (r) define the radial r de¬ 
pendence of the basis functions. The contribution of 
Ra,n,i,i',m' to the basis functions is much smaller than 
the contribution of R a ,n,i- The detailed shapes of the ra¬ 
dial functions R a ,n,i (0 and R a ,n,i,i',m' ( r ) are subject to 
the fitting algorithm. 

Step 1: First, electronic band structures e^ b (k) and 
wave functions are solved which serve as fitting tar¬ 
gets to the overall mapping algorithm 

H Ab (k) = sf b (k) (V’yk) ■ (2) 

The index j corresponds to the band index and k rep¬ 
resents a momentum vector in the first Brillouin zone. 
In principle, any method that is capable of solving band 
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updated basis 
functions and 
TB parameters 


FIG. 2: The process of mapping from ab-initio calculations to 
Tight Binding by which the ETB parameters and ETB basis 
functions are extracted iteratively. 


diagrams and explicit basis functions can provide these 
fitting targets. Throughout this work, however, hybrid 
functional calculations are performed for step lp2 

Step 2: In the second step, initial guesses for the ETB 
basis functions and ETB parameters are defined. During 
the fitting process, the ETB basis ‘Bmitiai is spanned by 
non-orthogonal functions {& a ,n,i,m ( r )} given by 

®a,n,L,m (r) = Yi >m (0, (j)) R a , n> i (r) . (3) 

The R a , n ,i (r) in Eq. © differ from the R a , n ,i (r) of the 
final basis functions in Eq. ([l]). The details of the ini¬ 
tial guesses for the diagonal and off-diagonal elements of 
the Hamiltonian H TB (k) are not essential for the over¬ 
all algorithm. Nevertheless, initial guesses that follow 
the framework of existing ETB parameter sets improve 
the overall fitting convergence. Urban et al. and Lu et 
al. discuss that interactions up to third nearest neigh¬ 
bors migh t be needed to exactly reproduce ab-initio re¬ 
sults.^® 1 'In contrast, we find that the interatomic inter¬ 
action elements of H TB (k) can be limited to first nearest 
neighbor interactions throughout this work while still re¬ 
producing ab-initio results very well. 

Step 3: The nonorthogonal basis functions $> a ,n,l,m (r) 
in position space are transformed into the Bloch repre¬ 
sent at ioiP^3> 0iTli 2 imik (r) 


© ortho, the wave functions 
formed into this basis 


V’/k ) °f step 1 must be trans- 


« P (k) IO = £ C j>a (k) |* a , k > , (5) 

Ot 


where 


c Jia (k) = (* aik |V$>, (6) 


C'j.a (k) ^^Q.k 


£ k ) -with the projection operator 
P( k) = £|* a , k ><* a , k [. (7) 


Equation ©> contains an approximation of the ab-initio 
wave functions in so far that the sum over a extends only 
over those orbitals that are included in the tight binding 
basis ©ortho- This basis and © or tho of similar ETB mod¬ 
els have much fewer basis vectors than the input ab-initio 
calculation. This rank reduction is a typical outcome of 
rectangular transformations such as P and is well known 
in the field of low rank approximations.® 5 

Step 4: Here, the quality of the ETB fitting is as¬ 
sessed. In this step, the band structures of the current 
ETB model eJ B (k) and the ab-initio input e Ab (k) are 
compared. If these sufficiently agree, the phases of the 
ETB wave functions are modulated to agree with the ab- 
initio ones and both wave functions are compared after 
that. The ETB Hamiltonian of step 2 is diagonalized in 
the basis ©ortho of step 3 to obtain ETB band structures 


sJ B (k) and eigen vectors 




H TB (k) (V’Jk ) = S J B (k) \ipj,u) > 


( 8 ) 


with 


|^Ik) = ££> (k) IdVk) - 


(9) 


To assess the quality of the ETB results is assessed, 
different fitness functions F e .F m and F^ are defined for 
energies, masses and wave functions respectively. The F e 
and F m are given by 


l^a.k) — ^a,n,Z,m,k (**) 

= ^ eX P I^ k ' ( R + T «)] (r - R - T a jj4) 

R 

where r a is the position of atom type a in the unit cell 
and the sum runs over all unit cells of the system with R, 
the position of the respective cell. To improve readabil¬ 
ity of all formulas in the Dirac notation, the indices of 
atom type and quantum numbers are merged into Greek 
indices a = (a, n, l, m). For the further steps, an orthog¬ 
onal basis ©ortho = {l^a.k)} is created out of the basis 
©initial with Lowdin’s symmetrical orthogonalization al- 
gorithm.^Since steps 4 and 5 are formulated in the basis 


F e = £Wj (k) I ej B (k) - ef b (k)f . (10) 

jM 


F m — £ ' 


m Ab — m TB 


->Ab 


(ii) 


where Wj (k) and iv rn are weights defined for each target. 

As a convention for wave functions phases, another set 
of ETB wave functions 

1tb\ 


4’Jk/ is introduced 
$j,k ) = £ V J,i (k) | £ k ) • 


(12) 
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The unitary transformation V (k) is defined by 

Vj,i (k) = 




A( 

k) 


(13) 


with 


x ^=JnT, 




V p £ 


(14) 


Here, the sum over p and q runs over all N ETB states 

with equivalent en- 


ipp .k and N ab-initio states 
,vrr; D o TB /k\ ~ ~Ab /l-\ flXo 


p.k J 

ergies e BB (k) 
equation holds 


s Ab (k) . With this transformation, the 




Ab 






= A(k), 


(15) 


for equivalent states. This phase adaption can only work 
if the ETB band structure is close enough to the ab-initio 
result. The ETB wave function fittness is given by 


f 4> = Y1 w t ( k ) |V$) ~ ) 

jM 


(16) 


The weights itif (k) are varying depending on respective 


$ 


,tb\ 


y,k J 


from 




- 4 £.) have in 


fitting focusses. Deviations of 

general two reasons: inadequate basis functions and/or 
eigenfunctions of a poorly approximated ETB Hamilto¬ 
nian. Therefore, can be estimated as 


*&)~ 






< 2 


+ 2 


i-PM |^> 


PM\^tt)- 






(17) 


The first right hand side term of the last equation de¬ 
scribes the deviation of the low-rank approximated ab- 
initio wave functions. This becomes obvious with the 
projector property P 2 (k) = P (k) 


i-PM | 


= m 


Ab 
'j, k 


i-PM 


^k) • 


(18) 

The second term on the right hand side of Eq. © con¬ 
tains information about the quality of the eigenfunctions 
of the approximate ETB Hamiltonian H TB (k). This is 
understandable when Eqs. ([5| and (12) are inserted into 
this term 


P{ k) 




ipi 


TB 


J-k/ 


2 — 2 Re 


E 


a,,4a m,i(kRa (k) 


(19) 


The fitness function F .represents the major improve¬ 
ment over the traditional ETB eigenvalue fitting (e.g. 
typically limited to energies and effective masses). All 
fitness functions are minimized by iterating over the steps 


3 and 4: the Slater Koster type parameters for the ETB 
Hamiltonian H TB (k) and the parameters of the radial 
ETB basis functions R a , n ,i (r) are adjusted for every it¬ 
eration of step3. 

Step 5: Once the fitness functions are small enough to 
cease the iterations, it is assumed that those eigenfunc¬ 
tions of the ETB Hamiltonian H TB (k) that were subject 
to the fitting are identical to the eigenfunctions of the ab- 
initio Hamiltonian H Ab (k) after a transformation A (k) 

l€k)«E^( k )|^)- (20) 

i 

This transformation A is determined by a singular value 
decomposition of the rectangular overlap matrix of ab- 
initio eigenstates with ETB eigenstates 

<4# l€k) = E ( k ) ( k ) w p.i ( k ) • ( 21 ) 

V 

The row index i runs over all ab-initio eigenstates - ex¬ 
ceeding those that served as fitting targets, whereas the 
column index j covers all the ETB eigenfunctions. The 
E and W are square and U is a rectangular matrix. The 
transformation A is then defined as 

A j,i (k) = E W 3,p ( k ) u l,i ( k ) ‘ (22) 

p 


A is constructed from relevant columns of a unitary trans¬ 
formation. Combining Eqs. (20) and § allows to deter¬ 
mine the Bloch periodic final basis functions 

I Kt) = E <i (k) A o,i (k) K£> • (23) 


The real space counterpart of 




final \ 
ct,k j 


is given by 


^fina.1 ( r _ R _ r ) 



III. RESULTS 

In this work, ab-initio level calculations of Si and GaAs 
systems were performed with VASPP^The HSE06 hybrid 
functional is used to produce reasonable band gaps in 
both the bulk and the UTB cases. In all HSE06 calcula¬ 
tions, a cutoff energy of 350eV is used. T-point centered 
Monkhorst Pack kspace grids are used for both bulk and 
UTB systems. The size of the kspace grid for bulk calcu¬ 
lations is a 6x6x6, while one for UTB is 6x6x 1. k-points 
with integration weights equal to zero are added to the 
original 6x6x6 or 6x6x1 grids in order to generate 
energy bands with higher k-space resolution. The spin 
orbit coupling is included in band structure calculations. 
Small hydrostatic strains up to 0.3% are introduced to 
adjust the bulk band gaps in order to match experimen¬ 
tal results. The lattice const used in this work is given 
by table |T| 
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FIG. 3: Band structure and density of states of bulk Si. ETB band structure agree with the HSE06 band structure (a), 
especially for bottom conduction bands (b) and top valence bands(c) around Fermi level. 



Bottom CB and top VB of bulk GaAs 



FIG. 4: Band structure and density of states of bulk GaAs. ETB band structure agree with the HSE06 band structure (a), 
especially for bottom conduction bands (b) and top valence bands(c) around Fermi level. 


A. Application to Bulk Materials 


For bulk Si and GaAs, fitting targets include the band 
structures of the lowest 16 bands (with spin degener¬ 
acy) along high symmetry directions, important effective 
masses and wave functions at high symmetry points such 
as r, L and X points. ETB basis functions in real space 
is reconstructed on 6 x 6 x 6 T center k space grid using 
Eq ([24). 


The band structures and DOS of bulk Si and GaAs 
(HSE06 vs ETB) are shown in Fig. [3] and [4] respectively. 
The band structures using existing Si and GaAs ETB 
parameter^ 1 27 are also shown in corresponding figures. 
The ETB band structures and DOS using parameters 
generated by this work show better agreement with the 


corresponding hybrid functional results compared with 
the existing parameterizations. For bulk Si, the existing 
parameterization shows a unexpected low s* band around 
5 eV above topmost valence bands. In the traditional fit¬ 
ting process, the s* band shows a strong preference for 
moving downward. 27 ^ Due to large number of parameters 
to be determined, traditional (energy-gap and effective- 
mass) based fitting procedures can find local minima in 
their fitness functions corresponding to wave functions 
significantly different from those predicted by ab-initio 
methods. The present method has the important advan¬ 
tage that optimization involves not only masses and gaps 
but also wavefunctions. Thus the ETB wavefunctions can 
be kept close to their ab-initio counterparts. For GaAs, 
the existing parameterization shows 2 eV higher s-type 
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Si 

GaAs 

a 0 

5.43 A 

a 0 5.6307A 



E s 

-2.803316 

E aa -8.063758 

E Sc 

-1.603222 

E p 

4.096984 

E Pa 3.126841 

Ep a 

4.745896 

E s . 

25.163115 

E a * 21.930865 

E** c 

23.630466 

E d 

12.568228 

E da 13.140998 

E dc 

14.807586 

A 

0.021926 

A a 0.194174 

A c 

0.036594 

Vsscr 

-2.066560 

V SaSc<T -1.798514 



Vs* s* CT 

-4.733506 

-4.112848 



V^s* a 

-1.703630 

V SaS * a -1.258382 

Vs^s* a 

-1.688128 

Vspcr 

3.144266 

V s „ Pca 3.116745 

Vs rp n CT 

2.776805 

Vs*pcr 

2.928749 

V s> ,.„ 1.635158 

Vs*p a CT 

3.381868 

V^dcr 

-2.131451 

V s „ dc a -0.396407 

Vs r dn. G 

-2.151852 

Vs* da 

-0.176671 

V <d ,.„ -0.145161 

Vs*d a cr 

-0.810997 

Vppcr 

4.122363 

V PaPoCT 4.034685 



Vpp7T 

-1.522175 

V PaP c7r -1.275446 



Vpdcr 

-1.127068 

V Va d a v -1.478036 

V Vr , d n <y 

-0.064809 

Vprf7T 

2.383978 

V Padc n 1.830852 

Vvrdn/K 

2.829426 

Vddcr 

-1.408578 

V dadc a -1.216390 



Vdci7r 

2.284472 

V dado n 2.042009 



Vdds 

-1.541821 

V dada s -1.829113 



H Si 

H c and H a 

Es h 

-3.056510 

E SHo 2.758428 

^ B H a 

-0.308397 

^SffS Si <T 

-4.859509 

V SH s a <r -2.960420 

S c 

a -3.151427 

ks/fPSi CT 

3.776178 

V SHoPa a 5.490764 

Gjjf a p c 

a 3.539284 

^SffS^CT 

0.0 

Vs H 0.0 

^s/f a sj 

CT -0.129904 

V SH dsi(T 

-0.007703 

V SHo d a a -1.727690 

a d c 

a -0.252733 

Ssi 

-0.276789 

5 a -0.266815 

5 C 

-0.586952 


TABLE I: Slater Koster type ETB parameters of bulk Si and 
GaAs, and passivation parameters of UTBs. 


low lying valence bands. The ETB parameters of bulk Si 
and GaAs are listed in table HI It can be seen from tables 
P and [Till the anisotropic hole masses by ETB show a 
remarkable agreement with HSE06 results. The principal 
authors of the previous works 11 27 explicitly pointed out 
that fitting hole masses had been very difficult with the 
previous methods. 

The orthogonal ETB basis functions Sfinai of Si, Ga 
and As atoms are shown in Fig. [bj The ETB basis func¬ 
tions are slightly environment dependent because they 
are orthogonal. Thus the ETB basis functions are not 
invariant under arbitrary rotations but invariant under 
symmetry operations within T d group, as pointed out by 
Slater and Koster PS It can be seen from Fig. [ 5 ] (a) to 
(f) that the s and p orbitals show s and p features near 
the atom. More complicated patterns in the area fur¬ 
ther away from the atom can be observed. These com¬ 
plicated patterns correspond to components with high 
angular momentums. The feature of orthogonal ETB 
basis function resembles the augmented basis functions 
used in ab-initio level calculations such as Augmented 
Plane Waves(APWs) and Muffin Tin Orbitals(MTOs). 
The orthogonal ETB basis functions have multiple angu¬ 
lar parts in each orbital as shown by Fig. [5](g),(h) and 
(i). The s, p and d type ETB basis functions are dom¬ 
inated by components with l = 0, 1 and 2 respectively. 
More than 90% for the s,p and d orbitals are comprised of 



Si 

targets 

TB Ref HSE06 

TB error (%) 

E g {T) 

3.399 

3.302 

3.244 

1.8 

E g (X) 

1.131 

1.142 

1.139 

0.2 

E g (L) 

2.383 

2.247 

2.188 

2.6 

A so 

0.047 

0.051 

0.052 

0.8 

nihhioo 

0.299 

0.281 

0.282 

0.097 

rrihh 110 

0.633 

0.566 

0.572 

0.977 

mhh.111 

0.796 

0.704 

0.714 

1.433 

mihioo 

0.232 

0.206 

0.204 

1.001 

mm 110 

0.165 

0.151 

0.149 

0.937 

mihin 

0.156 

0.143 

0.142 

0.927 

nisoioo 

0.266 

0.244 

0.242 

0.809 

nrisoiio 

0.266 

0.244 

0.242 

0.795 

msoiu 

0.267 

0.244 

0.242 

0.770 

rricxi 

0.887 

0.928 

0.857 

7.615 

nricxt 

0.225 

0.207 

0.215 

3.544 


TABLE II: Targets comparison of bulk Si. Critical band edges 
and effective masses at T, A' and L points by ETB and HSE06 
calculations are compared. 



GaAs 

targets 

TB Ref HSE06 

TB 

error(%) 

e 9 ( r) 

1.424 

1.418 

1.416 

0.2 

E g (X) 

1.900 

1.919 

1.910 

0.5 

E g (L ) 

1.707 

1.702 

1.708 

0.3 

A so 

0.326 

0.368 

0.367 

0.1 

nihhioo 

0.383 

0.310 

0.337 

8.510 

nihhuo 

0.667 

0.573 

0.619 

7.879 

nrihhiii 

0.853 

0.750 

0.813 

8.507 

niihioo 

0.085 

0.082 

0.083 

0.744 

mihuo 

0.078 

0.073 

0.074 

1.614 

niihin 

0.076 

0.071 

0.072 

1.715 

m S oioo 

0.166 

0.164 

0.160 

1.998 

ni S o 110 

0.166 

0.164 

0.160 

2.037 

msoiu 

0.166 

0.164 

0.160 

2.041 

nicioo 

0.068 

0.065 

0.067 

2.787 

nidio 

0.068 

0.066 

0.067 

2.790 

nidii 

0.068 

0.065 

0.067 

2.781 

m cX i 

1.526 

1.577 

1.480 

6.142 

rricxt 

0.177 

0.215 

0.204 

5.083 

men 

1.743 

1.626 

1.446 

11.055 

m C Lt 

0.099 

0.111 

0.136 

22.614 


TABLE III: Targets comparison of bulk GaAs. Critical band 
edges and effective masses at T, X and L from TB and HSE06 
calculations are compared. 


their l = 0, 1 and 2 components respectively. The excited 
s* type ETB basis functions have higher angular momen¬ 
tum and the / = 0 components have contributions of 60% 
to 70%. The second largest contribution in s* orbitals is 
the / component with l = 3. The / component attached 
to the s* orbitals have angular part equivalent to real 
space function xyz. This is a result of the existence of 
xyz-like crystal field near each atom in zincblende and 


























7 


Si - s 



(a) 


Ga - s 



Ga - py 


(b) 


(e) 


As - s 




(c) 




Angular quantum number I 



FIG. 5: Contours of selected ETB basis functions of Si((a),(d)), Ga ((b),(e)) and As ((c),(f)) atoms. (a),(b) and (c) correspond 
to the contours of s orbitals of Si, Ga and As atoms in x-y plane. (d),(e) and (f) correspond to the contours of p y orbitals of Si, 
Ga, and As atoms. (g),(h) and (i) show the contribution of different angular momentums in basis functions of Si ,Ga, and As 
atoms. The ETB basis functions of Si and GaAs are highly localized basis functions with one dominant angular momentum. 


diamond structures. 

B. Application to UTBs 

To validate the transferability of the ETB model, band 
structures and eigen functions of [001] UTBs passivated 
by Hydrogen atoms are calculated by both HSE06 and 
ETB models. The current calculations assume no strain 
in the UTBs. In the HSE06 calculations, charged hydro¬ 
gen atoms are used to passivate the dangling bonds of 
the surface atoms in GaAs UTBs. The surface As and 
Ga atoms are passivated by charged hydrogen atoms with 
3/4 ( denoted by H c ) and 5/4 ( denoted by H a ) elec¬ 
tron respectively. The charged hydrogen atoms neutral¬ 
ize most of the surface induced electric field in the UTBs. 
As a result, the charge distribution and local potential 
shows almost flat envelopes inside the UTBs. Small de¬ 
viation of potential can only be observed at the surface 
Si/Ga/As atoms. The nearly flat potential envelope sug¬ 
gests geometry dependent build-in potentials are needed 
only for surface atoms. Thus the comparisons between 


self-consistent hybrid functional calculations and single 
shot ETB calculations are fair. 

The HSE06 calculations show that the Hydrogen or¬ 
bitals contribute to the deep valence bands, thus Hydro¬ 
gen atoms are considered explicitly into the ETB Hamil¬ 
tonian of UTBs in this work. Is orbital is used as the 
ETB basis function for Hydrogen atoms. The explicit 
passivation model include extra Slater Koster type ETB 
parameters E ShS , V SHS m ^shp^^shs*u and Vsndcr- Fur¬ 
ther more, a geometry and element dependent potential 
6 is included for surface atoms. The onsite energies of 
the surface atoms are shifted by S. The onsite energy of 
the surface Ga atoms are thus E ao + Sh 0 ', and for sur¬ 
face As atoms, the onsite energies are E aa + 8 h c • Here 
the a stands for s 7 p,d and s* orbitals. ETB parameters 
of Si/GaAs in Si/GaAs UTBs are identical with the pa¬ 
rameters of unstrained bulk materials provided in section 

MM 

To determine the ETB parameters of H-passivation, 
band structures and real space wave functions of selected 





























(a) Si UTB (b) As terminated GaAs 
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FIG. 6: Band structures of 001 Si (a) and As terminated GaAs (b) UTBs by ETB agree with HSE06 band structures, 
demonstrating the bulk Si and GaAs ETB parameters are transferable to UTB cases. All UTBs contain 17 non-Hydrogen 
atomic layers(with thickness 4ao). 


(a) Si UTB (b) As terminated GaAs UTB 
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FIG. 7: Planar averaged real space probability amplitudes 
of lowest conduction and topmost valance states of 001 Si 
(a) and As terminated GaAs UTBs (b) by HSE06 and ETB 
calculations. With the real space TB basis functions, the re- 
alspace probability amplitudes of TB calculations show rea¬ 
sonable agreement with the HSE06 probability amplitudes. 
UTBs contain 17 non-Hydrogen atomic layers(with thickness 
4a 0 ). 


bands near the Fermi level of the UTBs are considered 
as fitting targets. The inclusion of wave functions as tar¬ 
gets serves the purpose of correcting possible problem¬ 
atic states. The target Si/GaAs UTBs contain 17 non- 
Hydrogen atomic layers. Parameters for Hydrogen atoms 
are also shown in table |T| In GaAs UTBs, As and Ga 
are passivated by Hydrogen atoms with different charge, 
thus the Hydrogen atoms have different onsite energies 
when different types of atoms are passivated. The Hydro¬ 


gen atoms bonding with As atoms are charged positively 
while the ones bonding with Ga atoms are charged neg¬ 
atively. Consequently, the H c which forms bond with 
As have a higher onsite energy than the H a which forms 
bond with Ga. 

Band structures of Si/GaAs UTBs are shown in Fig. [6] 
The ETB band structures match the HSE06 band struc¬ 
tures well for energies ranging from leV below the top¬ 
most valence bands to leV above the lowest conduction 
bands. Using the explicit ETB basis functions, ETB wave 
functions of UTBs with subatomic resolution are ob¬ 
tained and can be compared with corresponding HSE06 
wave functions. Planar averaged probability amplitudes 
of wave functions of the lowest conduction band and top 
most valence bands in Si/GaAs UTBs are shown in Fig. [7] 
It can be seen that not only the envelope but also de¬ 
tails in subatomic resolution of the ETB planar averaged 
\ijj| 2 show agreement with corresponding HSE06 results. 
On the other hand, Fig. [8] compares the ETB atom site 
resolved probability amplitudes among ETB models in 
present and previous works (Ref llll271) . The cations and 
anions in GaAs UTBs form different envelopes for all of 
the presented states. The lowest conduction and high¬ 
est valence states turn out to be well confined states 
in Si UTBs in all of the calculations. While, in GaAs 
UTBs, the lowest conduction states has significant con¬ 
tribution from the surface atoms. In Si ETB probabil¬ 
ity amplitudes by previous parametrizations show similar 
envelopes compared to the ETB and HSE06 probability 
amplitudes in this work. Fig. [8] (d) shows the problem¬ 
atic valence states in As terminated GaAs UTB by pa¬ 
rameters of previous work. The corresponding valence 
states by this work turn out to be a well confined ones. 
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FIG. 8: ETB atom site resolved probability amplitudes of Si 
((a),(c)), and As terminated GaAs ((b),(d)) UTBs using ETB 
parameters in this work and previous work. 111 ^® The ETB 
atom site probability using different parameters are qualita¬ 
tively similar in Si UTB, while the ETB atom site probability 
in As GaAs are more sensitive to the parameter sets and pas¬ 
sivation models, i.e. the valence states with parameters and 
passivation model by previous work are not confined. UTBs 
contain 17 atomic layers(thickness is 4do)- 


Envelopes of ETB valence states in GaAs UTB (As-t) 



FIG. 9: Comparison of ETB wave functions using different 
ETB parameters and passivation model, (a) and (b) use ETB 
parameters in ref. 1 ITl (c) and (d) use ETB parameters in this 
work, (a) and (c) correspond to implicit passivation model® 
(b) and (d) correspond to explicit passivation model. The 
ETB parameters with the explicit passivation model shows 
the most confined states, while the previous parameters and 
implicit passivation model lead to less confined states. 


To investigate this issue in more detail, in Fig. [9j ETB 
atom site resolved probability amplitudes for the top¬ 
most valence states of the four possible As-terminated 
GaAs UTBs are plotted: (a) previous parameter ^ 1 an d 
implicit passivatiorP^l (b) previous parameters and ex¬ 
plicit passivation; (c) new parameters and implicit pas¬ 
sivation; (d) new parameters and explicit passivation. It 
is clear that, for a given set of bulk parameters, the im¬ 
plicit passivation model leads to wavefunctions that are 
less-confined than those of the explicit passivation model. 
On the other hand, with the same passivation model, 
the ETB parameters by this works shows more confined 
top valence states than the existing ETB parameters. 
Thus the un-confined ETB state using the existing pa¬ 
rameter set and implicit passivation model appears to be 
due to both the bulk GaAs parameters and the passi¬ 
vation model. The implicit modeP s replaces the s- and 
p-orbitals of the surface atoms by sp3 hybrids and raises 
the energy of the dangling hybrids by S sp 3 = 30eV. The 
d- and s*-orbitals are left completely un-passivated, and 
the unconfined states of Fig. [ 9 ] (a) are only slightly af¬ 
fected by changing the value of S sp3 . The impact of al¬ 
ternate implicit passivation model to explicit passivation 
model is obvious by comparing sub-figures (a) to (b), as 
well as (c) to (d). To better understand the role of bulk 
parameters to this behavior, we experimented by reduc¬ 
ing the magnitude of the nearest-neighbor p a -d c coupling 
parameters in both sets as V Pa d a a —> V Pa d c a + 0.3eV, 
Vp a d 0 -K —• y Vp a d a TT ~ 0.3eU. Remarkably, in both cases 
the topmost valence-band state became much more con¬ 
fined. Bulk valence band wave functions in modified and 
original parameter sets tell the story: The general trend 
is that bulk sets which generate more p-like top of VB 
states give better confinement under passivation (and es¬ 
pecially implicit passivation) than do those with higher 
d-content. The reduction of |V), a d c0 -| and |Vp a d cir | lead to 
more p-like top VB states. Ga terminated case has less 
passivation problems because its top-of-VB bulk states 
have more contribution from the As atoms than from 
the Ga atoms. 


Fig. 10 shows the band gaps of the Si and GaAs [001] 
UTBs as functions of UTB thickness. With the ETB pa¬ 
rameters by this work, the ETB bandgaps of Si and GaAs 
UTBs with thickness from 0.5nm to 4nm agree well with 
the gaps by HSE06 calculations. The ETB bandgaps 
of Si UTBs using parameters from previous work also 
show good agreement with the HSE06 results. How¬ 
ever the ETB bandgaps of GaAs UTBs using parameters 
from previous work and implicit passivation model are 
of around 20% lower than the Hybrid functional results. 
The gaps of GaAs UTBs terminated with Ga and As 
atoms are very close in value for both Hybrid functional 
and ETB results in this work, however the gaps of GaAs 
UTBs terminated with Ga and As atoms by previous pa- 
rameterizations and implicit passivation model show 0.1 
to 0.2eV discrepancies. The band gap change in Si UTBs 
thicker than 3nm can be model by effective mass model( 
assuming parabolic E-k relation ). While in the GaAs 
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(a) Si [001] UTBs 



(b) As terminated [001] GaAs UTBs 



FIG. 10: Band gaps of Si UTBs (a) and As terminated UTBs (b) by HSE06 and ETB calculations. For the presented UTBs 
with thickness ranging from lnm to 4.5nm, the ETB band gaps have discrepancies of less than lOmeV compared with HSE06 
ones. The band gap changes by effective mass calculation show agreement with HSE06 for Si UTBs thicker than 3nm. While the 
effective mass calculations has obvious discrepancies for all GaAs UTBs. The HSE06 and ETB calculations using parameters 
by this work consider Hydrogen atoms explicitly, while the ETB calculations using parameters by previous work is based on 
implicit passivation model. 12 


UTBs, the discrepancies between effective mass calcula¬ 
tions and HSE06 or TB calculations are obvious for all 
GaAs UTBs presented, suggesting the non-parabolic fea¬ 
ture of the GaAs valleys have significant impact to GaAs 
nano structures. The gaps by previous parameterization 
with implicit passivation model of As terminated GaAs 
UTBs has lower confined energies due to the unconfined 
valence states. 

IV. CONCLUSION 

It has been shown that the existing ETB parameteriza¬ 
tion together with the implicit passivation model gives 
unphysical states in As terminated GaAs UTB calcula¬ 
tions. A more reliable technique of ab-initio mapping 
which generates ETB parameters and basis functions 
from ab-initio is developed. The ab-initio mapping pro¬ 
cess is applied to both bulk Si and GaAs. Slater Roster 
type ETB parameters within 1st nearest neighbour ap¬ 
proximation and highly localized ETB basis functions are 


obtained. The ETB parameters and basis functions of 
Si and GaAs are validated in corresponding UTB sys¬ 
tems with passivation models that consider Hydrogen 
atom explicitly. Band gaps in Si and GaAs UTBs with 
different thickness are also calculated by HSE06, ETB 
and effective mass model. Compared with the existing 
ETB parameterizations and implicit passivation model, 
the ETB calculations in this work show good agreements 
with HSE06 calculations in both band structures and 
wave functions. This work shows that the ETB parame¬ 
ters by ab-initio mapping have good transferability. The 
mapping method developed here significantly reduces the 
uncertainty in both bulk and passivation models. 
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